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ABSTRACT 

The clustering of galaxies relative to the underlying mass distribution declines with cosmic 
time for three reasons. First, nonlinear peaks become less rare events as the density field evolves. 
Second, the densest regions stop forming new galaxies because their gas becomes too hot to 
cool and collapse. Third, after galaxies form, they are subject to the same gravitational forces 
as the dark matter, and thus they tend to trace the dark matter distibution more closely with 
time; in this sense, they are gravitationally "debiased." In order to illustrate these effects, we 
perform a large-scale hydrodynamic cosmological simulation of a ACDM model with f^o = 0.37 
and examine the statistics of <5*(r, z), the density field of recently formed galaxies at position r 
and redshift z. We find that the bias of recently formed galaxies 6* = (S 2 ) 1 ^ 2 / (S 2 ) 1 ^ 2 , where 
S is the mass overdensity, evolves from 6* ~ 4.5 at z = 3 to ~ 1 at z = 0, on 8 h^ 1 Mpc 
comoving scales. The correlation coefficient = (55*) / (S 2 ) 1 ^ 2 (S 2 ) 1 ^ 2 evolves from r* ~ 0.9 at 
z = 3 to r* ~ 0.25 at z = 0. That is, as gas in the universe heats up and prevents star formation, 
the star-forming galaxies become poorer tracers of the mass density field. We show that the 
linear continuity equation is a good approximation for describing the gravitational debiasing, 
even on nonlinear scales. The most interesting observational consequence of the simulations 
is that the linear regression of the galaxy formation density field on the galaxy density field, 
6» g r* s — (5 r 5 g )/{5 2 ), evolves from about 0.9 at z = 1 to 0.35 at z = 0. Measuring this evolution, 
which should be possible using the Sloan Digital Sky Survey, would place constraints on models 
for galaxy formation. In addition, we evaluate the effects of the evolution of galaxy formation 
on estimates of f2 from cluster mass-to-light ratios, finding that while Q(z) increases with z, the 
estimated Q cs t(z) actually decreases. This effect is due to the combination of galaxy bias and 
the relative fading of cluster galaxies with respect to field galaxies. Finally, these effects provide 
a possible explanation for the Butcher-Oemler effect, the excess of blue galaxies in clusters at 
redshift z ~ 0.5. 
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Motivation 



Recent observations of the evolution of galaxy clustering have explored both the regime z < 1, and 
have probed the distribution of galaxies at z 3. The low-redshift works include a wide- field (4° x 4°) 



/-selected angular sample ( Postman et al. 1998 ), Canada- France Redshift Survey galaxies with redshifts 
between < z < 1.3 ( [Lc Fevre et al. 1996 ), the Hubble Deep Field ( |Connolly, Szalay, fc Brunner 1998j ), 
field galaxies from the Canadian Network for Observational Cosmology (CNOC) survey ( Shepherd et 



. 1997), the Hawaii if-selected redshift survey (Carlberg et al. 1997), and a redshift survey performed on 



the Palomar 200-inch telescope (Small et al. 1999). The results of these surveys are still rather discrepant; 
they don't yet give a consistent picture for how galaxy clustering evolves with time. The evolution of the 
relative clustering of old and and young galaxies is also of interest, because it yields direct information on 



where galaxy formation occurs as a function of redshift. Lc Fevre et al. (1996) made a simple attempt at 
measuring this evolution by separating their sample by color. They found that red and blue galaxies have 
comparable correlation amplitudes at z > 0.5, although red galaxies are more clustered today, by a factor 
of about 1.5 (in agreement with other local estimates, such as Davis fe Geller 1976] , |Giovanelli, Haynes, fc 
Chincarini 1986, Santiago & Strauss 1992, Loveday et al. 1995, Loveday et al. 1996, Hermit et al. 1996, and 



Guzzo et al. 1997). The evolution of the density-morphology relation in clusters also contains information 
about where galaxies form as a function of redshift, showing that at z ~ 0.5 there is significantly more star 



formation in clusters than there is today (Butcher & Oemler 1978, 1984; Couch fc Sharpies 1987; Dressier 



fc Gunn 1992j ; |Couch et al. 1994[ |Oemler, Dressier, fc Butcher 1997| ; |Dresslcr et al. 1997[ |Poggianti et al 



1999| ). 



Measurements of the clustering of Lyman-break objects (LBOs) at z ~ 3 ( [Btcidcl et al. 1998 ) give a 
more unambiguous indication of the evolution in clustering. In particular, measurements of the amplitude 



of counts- in-cells fluctuations ( Adelberger et al. 1998 ) or the angular autocorrelation function ( Giavaliscc 
et al. 1998) suggest that galaxies were as strongly clustered in comoving coordinates at z ~ 3 as they 



are today. If the LBOs are unbiased tracers of the mass density field, these results contradict the widely 
accepted gravitational instability (GI) model for the formation of large-scale structure, unless one assumes 
an unacceptably low value of Oo in order to prevent the growth of the clustering of galaxies between z = 3 
and today. Therefore, the objects observed at high redshift are probably more highly "biased" tracers of the 
mass density field than are galaxies today. That is, they have a high value of b g = cr g /a, where a g = {Sg) 1 ^ 2 
is the rms galaxy density fluctuation and a = (S 2 ) 1 ^ 2 is the rms mass fluctuation. The counts-in-cells 
analysis of Adelberger et al. (1998) suggests that at z = 3 the bias of LBOs is b g ~ 2bg (for f2o = 0.2, 
= 0), b g ~ Abo (for f^o = 0.3, ft a = 0.7), or b g ~ 6&o (for Slo = l)i where bo is the bias of galaxies today. 

How does the observed clustering since z ~ 3 relate to the underlying mass density fluctuations? From a 
theoretical perspective, the bias decreases with time because of three effects. First, at early times, collapsed 
objects are likely to be in the highest peaks of the density field, since one needs a dense enough clump 
of baryons in order to start forming stars. Such high-er peaks are highly biased tracers of the underlying 
density field ( poroshkevich 1970| ; [Kaiser 1984] ; [Bardecn et al. 1986j ; |Bond et al. 1991| ; |Mo fc White 1996| ). 



As time progresses, lower peaks in the density field begin to form galaxies, which are less biased tracers of 
the mass density field. Second, the very densest regions soon become filled with shock-heated, virialized 
gas which does not easily cool and collapse to form galaxies ( Blanton et al. 1999| ) . These two effects cause a 
shift of galaxy formation to lower density regions of the universe. This shift is evident in the real universe: 
young galaxies are not found in clusters at z = 0. To keep track of these effects, we define the bias fr* and 
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correlation coefficient of galaxy formation as 



1/2 



{<r- 



1/2 



and 



(SJ) 



l/2,, 2 >l/2' 



(1) 



where S = p/(p) — 1 is the overdensity of mass and J* = p*/{p*} — 1 is the overdensity field of "recently 
formed" galaxies, each defined at some smoothing scale. By "recently formed," we mean material that has 
collapsed and formed stars within the last 0.5 Gyrs. Alternatively, one can consider 6* as the overdensity 
field of galaxies weighted by the star formation rate in each — for the purposes of this paper we consider 
the formation of a galaxy equivalent to the formation of stars within it. As functions of redshift, b*(z) and 
T*(z) are essentially differential quantities which keep track of where stars are forming at each epoch. They 
both decrease with time for the reasons stated above. Predictions of b*(z) and r*(z) are complementary 



to predictions of the star formation rate as a function of redshift (Madau et al. 1996; Nagaminc, Cen fe 



Ostrikcr 1999 ; Somerville fc Primack 1998 ; Baugh et al. 199E ); whereas those studies examine when galaxies 



form, we examine where they form. 

The third effect is that after galaxies form, they are subject to the same gravitational physics as the 
dark matter; thus, the two distributions become more and more alike, debiasing the galaxies gravitationally 
(Fry 1996; Tegmark & Peebles 199? ; see Section 0). The history of and r*(z) convolved with effects of 



gravitational debiasing determines the history of the bias of all galaxies, quantified as 



m 1 ' 2 



,1/2 



and 



l/2/X2\V2 ; 



(2) 



where 8 g = p g /(p g ) — 1 is the overdensity of galaxies. Here, b g {z) and r g (z) are integral quantities, which 
refer to the entire population of galaxies. Until recently, most papers on this subject have considered only 
the effect of bias (b g ) on large-scale structure statistics. However, observational and theoretical arguments 



suggest that scatter in the relationship between galaxies and mass may also be important (Tegmark & 



Bromley 1998; Blanton et al. 1999); treating the relationship as purely deterministic when analyzing 



peculiar velocity surveys or redshift-space distortions of the power spectrum can produce observational 
inconsistencies ( Pen 199"g| ; Dckel fc Lahav 1998 ). Thus, even at lowest order, these treatments need to 
include consideration of the correlation coefficient r g . 

Finally, we will also define the cross-correlation of all galaxies and of young galaxies as 



b =m 



2\l/2 



^ " <0 1/2 



and 



<«> 1/a <*2> 1/a ' 



(3) 



These quantities are important because they are potentially observable, by comparing number-weighted 
and star-formation weighted galaxy density fields. 

Previous theoretical efforts have focused on the evolution of the clustering of all galaxies, or equivalently, 
the evolution of the integral quantities b g {z) and r g (z). There has been notable success in this vein in 
explaining the nature of the LBOs. In the peaks-biasing formalism, the halo mass of collapsed objects 



determines both their number density and their clustering strength ( Doroshkevich 197C; Kaiser 1984 



Bardeen et al. 1986[ [Bond et al. 1991 ; [Mo fc White 1996; ); interestingly, halos with masses > 10 12 Af G 



at z = 3 have both a number density and clustering strength similar to those of LBOs, in reasonable 



cosmologies (Adelberger et al. 1998), a result verified by A-body simulations (Wechsler et al. 1998; 
Kravtsov & Klypin 1998). This fact is reflected in the results for high-redshift clustering of semi- analytic 
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galaxy formation models, which also find that galaxies at z ~ 3 are highly biased (Kauffmann, Nusser, & 



Stcinmetz 1997: Kauffmann et al. 1999; Somerville, Primack & Faber 1998 



Baugh et al. 1998). Furthermore, 



hydrodynamic simulations using simple criteria for galaxy formation (Evrard, Summers, & Davis 1994; 
Katz, Hcrnquist, fc Weinberg 1998 ; |Cen fc Ostrikcr 1998a ) find that the distribution of galaxies at z = 3 



is indeed biased with respect to the mass distribution, and that the clustering strength of galaxies depends 
only weakly on redshift. 

While the issue of b g (z) and r g (z) is the one which observations are currently best-suited to address, 
perhaps other statistics can constrain the nature of galaxy formation more powerfully. After all, the 
fundamental prediction of theories for galaxy formation, such as ours and such as the semi-analytic models 
mentioned above, is the location of star formation as a function of time. Important information about 
galaxy formation may be lost if one examines only the integral quantities, which have been affected by the 
entire history of galaxy formation convolved with their subsequent gravitational evolution. Therefore, we 
focus here on the evolution of the large-scale clustering of galaxy formation, defined as the formation of the 
galaxy's constituent stars — that is, the evolution of 6* and r*. Now is an opportune time to investigate 
such questions, since large redshift and angular surveys which can examine this and related properties of 
galaxies are imminent. 

In this paper, we measure the clustering of galaxy formation in the hydrodynamical simulations of 
Cen fc Ostrikcr (1998a^ , which we describe in Section ^. In Section ||, we consider the properties of galaxy 
formation in these simulations as a function of time, studying the first two effects discussed above and 
the redshift dependence of 6* and r* . In Section ||, we discuss the gravitational debiasing, the third effect 
mentioned above. In Section ^, we discuss observable effects of the trends in galaxy formation found in the 
simulations. Especially important is the strong evolution of the cross-correlation between the star-formation 
weighted density field and the galaxy density field; this evolution should be observable with surveys such as 
the Sloan Digital Sky Survey (SDSS; Gunn fc Weinberg 1995| ). In addition, the evolution of the location of 
galaxy formation in the universe may be related to the observed Butcher-Oemler effect and may affect the 
mass-to-light ratios of clusters. We conclude in Section ||. 



Simulations 



The work of Ostriker fc Steinhardt (1995)| motivated the choice of a flat cold dark matter cosmology for 



the simulation used in this paper, with £!q = 0.37, tt\ = 0.63, and fib = 0.049. Recent observations of high 



redshift supernovae have lent support to the picture of a flat, low-density universe (Pcrlmuttcr et al. 1997 



Garnavich et al. 1998). Great uncertainty remains, however, and future work along the lines of this paper 
will need to address different cosmologies. The Hubble constant was set to Hq = 100 h km s _1 Mpc -1 , 
with h = 0.7. The primordial perturbations were adiabatic and random phase, with a power spectrum 
slope of n = 0.95 and amplitudes such that <7g = 0.8 for the dark matter at z = 0, at which time the 
age of the universe is 12.7 Gyrs. We use a periodic box 100 h^ 1 Mpc on a side, with 512 3 grid cells and 
256 3 dark matter particles. Thus, the dark matter mass resolution is about 5 x 10 9 h^ 1 Mq and the grid 
cell size is ~ 0.2 hr 1 Mpc. The smallest smoothing length we consider is a 1 h~ x Mpc radius top hat, 
which is considerably larger than a cell size. On these scales and larger, the relevant gravitational and 
hydrodynamical physics are correctly handled. On the other hand, subgrid effects such as the fine grain 
structure of the gas and star formation may influence large-scale properties of the galaxy distribution. As 
we describe here, we handle these effects using crude, but plausible, rules. 
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Cen & Ostriker (1998a) describe the hydrodynamic code in detail; it is similar to but greatly improved 
over that of Cen & Ostriker (1992a, b). The simulations are Eulerian on a Cartesian grid and use the 
Total Variation Diminishing method with a shock-capturing scheme (Jameson 1989). In addition, the code 
accounts for cooling processes and incorporates a heuristic galaxy formation criterion, whose essence is as 
follows: if a cell's density is high enough, if the cooling time of the gas in it is shorter than its dynamical 
time, if it contains greater than the Jeans mass, and if the flow around that cell is converging, it will have 
stars forming inside of it. The code turns the baryonic fluid component into collisionless stellar particles 
(hereafter "galaxy particles") at a rate proportional to mb/tdyn, where rrib is the mass of gas in the cells 
and £dyn is the local dynamical time. These galaxy particles subsequently contribute to metal production 



and the background ionizing UV radiation. This algorithm is essentially the same as that used by Katz 



Hcrnquist, fc Weinberg (1992)| and Gncdin (1996a, b). The masses of these galaxy particles range from 



about 10 to 10 Mq. Thus, many galaxy particles are contained in what would correspond to a single 
luminous galaxy in the real universe. However, rather than grouping the particles into galaxies, in this 
paper we simply define a galaxy mass density field from the distribution of galaxy particles themselves. 
More details concerning the galaxy formation criteria and their consequences are given in Appendix [X| 

We will examine the results of the simulations at four output times: z = 3, z = 1, z = 0.5 and 
z — 0. At each output time, we will consider the properties both of all the galaxy particles and of those 
galaxy particles formed within the previous 0.5 Gyrs. We will take these "recently" formed particles to be 
representative of the properties of galaxy formation at that output time, and label their overdensity field 
5* = p*/(p*) — 1, the "galaxy formation density field." 



3. Time dependence of Galaxy Formation 

In this section, we study the relationship between the mass density, the gas temperature, and the 
galaxy formation density fields at each epoch. We will show the evolution of the the quantities 6* and r* 
with redshift, and demonstrate that this evolution is largely due to the changing relationship between mass 
density and temperature. 

Let us begin by considering the evolution of the temperature and mass density fields as a function of 
redshift. At each of our four output times, we smooth the fields with a top hat of 1 /i -1 Mpc comoving 
radius. In Figure [l] we show the joint distribution of 1 + 8, the normalized mass density, and T, the gas 
temperature, as a function of redshift. We could instead have considered the distribution of 1 + 8b, the 
baryonic density; the results would have been nearly identical, because the dark matter and baryonic 
overdensities are closely related on these scales. Most of the gas, that is, does not cool, fall to the central 



parts of galaxies, and form stars, but instead follows the same distribution as the dark matter (Cen & 



Ostriker 1998b] ). As the universe evolves, the variance in density increases. In addition, as gas falls into the 



potential wells, it shocks and heats up to the virial temperature necessary to support itself in the potential 
well, raising the average temperature of gas in the universe. Thus, volume elements above the mean density 
move up and to the right on this diagram. These results are consistent with previous work with the same 



code (for instance, Cen, Gnedin, & Ostriker 1993) and with other codes (for a summary, see Kang ei 
al. 1994 ). Note that on 1 h" 1 Mpc scales, the importance of cooling physics to the gas temperature is not 
immediately apparent; in Appendix A, we show this same figure at the scale of a single cell, which shows 
quite clearly the effects of including the radiative physics of the gas. 

We now consider the equivalent plot for the galaxy formation density field. Figure @ shows the 
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conditional mean of the galaxy formation density field 1 + <5» as a function of 1 + 6 and T: (1 + 5* |1 + 8, T). 
At any fixed density, the star formation rate declines as temperature rises. It is clear that the densest 
regions cannot cool and collapse at late times because they are too hot. To illustrate what Figure || 
implies about the relationship between galaxy formation and mass density, Figure [| shows the conditional 
probability P(l + 5*|1 + 6) as the logarithmic grey scale and the conditional mean (1 + <5»|1 + S) as the 
solid line. The dashed lines arc la limits around the mean. At high redshift, the high density regions are 
still cool enough to have clumps of gas in them cooling quickly; in this sense, the temperature on these 
scales is not particularly important, making the relationship between galaxy formation and mass density 
close to deterministic. At low redshift, many regions become too hot to form new galaxies; in particular, 
galaxy formation is completely quenched in the densest regions, causing a turnover in (1 + 5*|1 + 6). In 
addition, at low redshifts, the galaxy formation rate depends strongly on temperature, adding scatter to 



the relationship between galaxy formation and mass (Blanton et al. 1999) 



The most obvious consequence of this evolution is that the bias of galaxy formation is a strong function 
of time. Figure ^ shows 6* and r* (defined in Section |l|) as a function of redshift for top hat smoothing 
radii of R — 1 and R — 8 h^ 1 Mpc. The bias of galaxy formation 6* is clearly a strong function of redshift, 
especially at small scales. The decrease in the bias is due to the fact that the galaxy formation has moved 
to lower a peaks and out of the hottest (and densest) regions of the universe, as described in Section [l| 

We can compare these results to those of the peaks-bias formalism for the bias of dark matter haloes. 
If we assume that the distribution of galaxy formation is traced by M > 10 12 Mq halos (that is, halos 
with typical bright galaxy masses) which have just collapsed, then according to the peaks-bias formalism 



(Doroshkcvich 1970; Kaiser 1984; Bardcen et al. 1986; Bond et al. 1991; Mo fc White 1996), the bias of 



mass in such halos with respect to the general mass density field is: 

where S c ~ 1.7 is the linear overdensity of a sphere which collapses exactly at the redshift of observation, 
and cr(M, z) is the rms mass fluctuation on scales corresponding to mass M, again at the redshift of the 
observations. Note that this formula differs from the standard one for the bias of all halos which have 
collapsed before redshift z. The dotted line in the upper panel of Figure |J shows the resulting bias. Thus, 
the prediction for 6* of the pure halo bias model is not much different from that of our hydrodynamic 
model on large scales, even though the former does not include all the physical effects of the latter. That 
is, measuring 6*(z) alone cannot easily distinguish the two models. 

More interesting is the evolution of the correlation coefficient r* as a function of redshift, shown in 
the lower panel of Figure ||. At all scales, the correlation coefficient falls precipitously for z < 1, when the 



densest regions become too hot to form new galaxies. On the other hand, in the case of peaks-biasing, Mo 



fc White (1996) find a scaling between the halo and mass cross- and autocorrelations that implies ~ 1 



on large scales at all redshifts^- This is an important difference between the predictions of the models of 



Mo & White (1996) and ours, which is presumably due to the more realistic physics used in our model. 
Semi-analytic models can predict the evolution of r* ; it will be interesting to see how dependent this statistic 
is to the details of the galaxy formation model. We discuss a way to get a handle on observationally in 
Section ||. 



'On scale s of 1 h 1 Mpc, things a re more complicated, because 10 12 Mq halos so close together actually overlap in Lagrangian 
space — see jheth fc Lemson 1998[ 
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Gravitational Evolution of Bias 



After galaxies form, they fall into potential wells under the influence of gravity. Because the acceleration 
on galaxies is the same as that on the dark matter, this gravitational evolution after formation will tend 



to bring both the bias b g and the correlation coefficient r g closer to unity, as described by Fry (1996 



and Tcgmark fc Peebles (1998) . The evolution of the bias, then, is determined by the properties of 
forming galaxies (outlined in the previous section) and how those properties evolve after formation. 
Bere we investigate this process and show that the linear approximations of Fry (1996)| and Tcgmark & 



Peebles (1998)| describe this evolution well, even in the nonlinear regime. 



4.1. Evolution of the Clustering of Coeval Galaxies 

I ~~l 

Fry (1996) and Tegmark & Peebles (1998) rely on the continuity equation to study the evolution of the 

density field of galaxies formed at a given epoch. In this subsection, subscript c will refer to this coeval set 
of galaxies. Both galaxies and mass satisfy the same equation, which is to first order 

<j + V-v = 0, (5) 

assuming that galaxies are neither created nor destroyed. Since we are dealing with a coeval set of galaxies, 
they cannot be created by definition. Although it is possible to disrupt or merge galaxies in dense regions 
( Kravtsov fc Klypin 1998| ), we are following the stellar mass density in this paper, not the galaxy number 



density, and thus do not consider these effects. 

Assuming that there is no velocity bias between the volume-weighted velocity fields of mass and of 
galaxies (which is borne out by the simulations), it follows that S = S c . This statement means that on 
linearly scaled plots of S c against 5, volume elements are constrained to evolve along 45° lines. We test this 
prediction on 30 hr 1 Mpc top hat scales, which should be close to linear, in Figure ^ Bere we consider 
a burst of galaxies formed at z = 1. The solid lines are (S C \S) at z = 1, z = 0.5 and z = 0. The dashed 
lines are the predictions extrapolated from z — 1 as above, to z = 0.5 and z = 0, respectively. Notice use 
the growth factor of mass, a(z = 0)/a(z — 1), which on these scales is close to its linear theory value, to 
determine AS, and assume that AS C is the same. Evidently the linear continuity equation works on these 
scales. We also test this prediction on the nonlinear scale of 1 h^ 1 Mpc in Figure [| We find that it works 
extremely well in the range 3 < S < 100; in higher and lower density regions, nonlinear corrections to 
Equation (^) are clearly important. 
i i 

Tegmark & Peebles (1998) use the continuity equation to make further predictions about the bias 

properties of coeval galaxies. In particular, given a "bias at birth" of b c (zo) and a "correlation coefficient at 
birth" of r c (zo), one can express b c (z) and r c (z) in terms of the linear growth factor relative to the epoch of 
birth D(z)/D : 

M*)r„(*) = 1+ b ^)r c (z )-l 
D(z)/D a 

,y , _ [(1 - D(z)/Dq) 2 - 2(1 - D(z)/D )b c (zo)r c (z ) + b 2 c (z )} 
c[Z) " (D(z)/D r [) 

In Figure [?] we show as the solid lines the values of b c and r c for the burst of galaxies at z = 1, at each output 
time. The dashed lines are the predictions based on z = 1 for the results at z — 0.5 and z — 0. Notice that 
on small scales, the prediction for b c remains extremely good. For r c , the linear prediction underestimates 
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the degree to which galaxies and mass become correlated on small scales; essentially, nonlinearities enhance 
the rate at which initial conditions are forgotten. 

We can compare the effects of debiasing to the decline in the correlation between galaxy formation and 
mass with redshift. In Figure || we show 6,(2) and r*(z) at 8 hr 1 Mpc scales, from Figure |[ as the solid 
lines. The dotted lines originating at each output redshift show how the galaxies formed at that redshift 
would evolve according the continuity equation. Naturally, b c [z) declines while r c {z) rises. Note that &*(z) 
falls more quickly than debiasing can occur, which is why old galaxies at redshift zero are still more highly 



biased than young galaxies (Blanton et aL 1999), despite having had more time to debias 



4.2. Evolution of the Clustering of the Full Galaxy Population 

These results show that the decrease in the bias of the galaxy density field with time found by Cn\ 



fc Ostrikcr (1998a)| must be partly due to gravitational debiasing (described in this section) as well as the 
decrease in 6* (described in Section [|) . We investigate this process here by showing the evolution of the 
bias of all the galaxies in Figure ^ as the solid lines, for a tophat smoothing of 8 hr 1 Mpc. The evolution of 
the bias of galaxy formation, from Figure^, is shown as the dashed lines. Notice that b g (z) and 6*(z) are 
nearly the same, even though r g {z) and r*(z) differ considerably; r*(z) plummets after 2 = 1, while r g (z) 
remains close to unity. About half the star formation in the simulation occurs at z < 1. Although this 
recent star formation has a low r„, and is thus poorly correlated with the mass distribution, this has little 
effect on r g , partly due to the gravitational debiasing. 



Using the approximations of Tegmark fc Peebles (1998) , we can reconstruct b g (z) and r g (z) from the 



properties of galaxy formation as a function of redshift. We get equations similar to those of equation 
now integrated over the star formation history: 

b g (z)r g (z) = 1 + - * , / dz' [K(z')r^z') - 1] D{z') dG 



D{z)G{z) J l * v j v 1 dz r 



b%z) = bl(z)rl(z) + ^ (7) 



where D(z) is the linear growth factor, G(z) is the cumulative mass of stars formed by redshift z, and a 
is the rms mass fluctuation. Calculating b g (z) requires first evaluating b g {z)r g {z) substituting it into the 
expression given. In the second term, the quantity s 2 is defined as: 



,2 



S 



1 f f ,„,,dG 



G 2 {z) J dz 1 J dz 



hi 



dz'— / dG" — (5 ± ( z i)6 ± (z")), (8) 



where <5j_ = i5» — 6*r»<5 is the component of the galaxy formation field which is uncorrelated with the 
local density.^ These equations are simply a result of applying the continuity equation (Equation ||) to 
the case of continuous star formation, rather than an instantaneous burst, as in Equation (^|). Figure || 
shows bg(z), r g (z), and b g (z)r g (z) for all galaxies reconstructed from the properties of our galaxy formation 
fields. b g {z)r g {z) represents the linear regression of the galaxy density on the mass density. Note that the 
reconstruction of this quantity is quite accurate; on the other hand, the reconstructions of b g (z) and r g (z) 
separately, which require knowledge of (6±(z')5±(z")}, have significant errors. These errors are mostly due 



Note that the definition of fc* in Tegmark & Peebles (1998) is equivalent to our quantity b„r 
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to the poor time resolution of our output j^| The continuity equation will prove useful to future theoretical 
and observational work, as we discuss in Section pi 



5. Observing the Evolution of Galaxy Formation 

Here we discuss the observational consequences of the evolution of the spatial distribution of galaxy 
formation. First, we consider observable properties of the star formation density field. Second, we consider 
the resulting properties of galaxy clusters. 



5.1. Star Formation Density Field 

We noticed in the last section that the correlation coefficient between star forming galaxies and mass 
should decrease considerably with time. We cannot currently observe this decrease directly in the real 
universe]^ However, from Figure we note that the galaxy distribution as a whole does correlate well with 
the mass, in the sense that r g (z) ~ 1 at all redshifts. Therefore, we should be able to detect the evolution 
of r*(z) by cross-correlating galaxy formation with the distribution of all galaxies. In Figure ^ we plot 
b^ g = <7*l<7 g and r* g = {5*6 g ) / o~*o~ g between the galaxy formation density field and all galaxies as a function 
of redshift, at 8 /i -1 Mpc scales. While the evolution in 6* g is rather weak, the evolution in r* 9 is striking, 
and should be observable. 

Observationally, measuring this correlation will require mapping the density field of star formation in 
the universe, which has not yet been done. Given that the fundamental prediction of galaxy formation 
models is the location of star formation as a function of time, such a map would be extremely useful. The 
spectral coverage of the Las Campanas Redshift Survey (LCRS) includes a star formation indicator ([Oil] 



equivalent widths; Schectman et al. 1996; Hashimoto et al. 1998); thus, one can measure this correlation 



at low redshift using this data. Using the LCRS, Tegmark & Bromley (1998) have addressed a related 



question by measuring the correlation between the distributions of early and late spectral types (as classified 



by Bromley et al. 1998|); they found r ~ 0.4-0.7, in qualitative agreement with our findings here, if one 



makes a correspondence between spectral type and star formation rate. They can calculate r because 
they understand the properties of their errors, which are due to Poisson statistics; however, in the case of 
measuring spectral lines to determine the star formation rate in a galaxy, the errors are much greater and 
more poorly understood. Therefore, since 5 g is likely to be much better determined than 5*, the quantity 
of interest will not be r* g , but the combination b* g r* g — (S*5 g )/o~ g , which is the linear regression of star 
formation density on galaxy density and is independent of cr*. Future redshift surveys such as the SDSS will 
probe redshifts as high as z ~ 0.2 — 0.3 and will have spectral coverage which will include measures of star 
formation such as [Oil] and Ha. Figure || indicates that 6* g r^ g ~ 0.35 at z = and b* g r„ g ~ 0.5 at z = 0.25 
(linearly interpolating between z = and z = 0.5), a difference that should be measureable. A cautionary 
note is that the linear interpolation in this range may not be accurate, and that simulations with finer time 



3 The results of these simulations take large amounts of disk space, and thus when they were run we did not attempt to store 
more than a handful of output times. It would be impractical to rerun the simulations again for the sole purpose of obtaining 
more output redshifts. 



4 Weak lensing studies (Mellier 1999) and peculiar velocity surveys (Strauss & Willick 1995) do probe the mass distribution 



and may prove useful for studying the evolution of b g (z) and r g {z) directly. 
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resolution are necessary to make more solid predictions. 

In addition, the SDSS photometric survey will probe to even higher redshift. Using the multiband 
photometry, one will be able to produce estimates both of redshift and of spectral type. These measurements 
will allow a calculation of the angular cross-correlation of galaxies of different types or colors as a function 
of redshift out to z ~ 0.5 — 1; Figure |l0| shows that this evolution should be strong. The drawback to this 
approach is that one cannot make a quantitative comparison between our models and the cross-correlations 
of galaxies of different colors, because we can not yet reliably identify individual galaxies in the simulations 
and because there is not a one-to-one correspondence between galaxy color and ages. Needless to say, the 
observational question is still worth asking, because quantitative comparisons will be possible in the future 
using hydrodynamical simulations at higher resolution, and are even possible today using semi-analytic 
models (Kauffmann et al. 1999; [Bomcrvillc, Primack fc Fabcr 1998; Baugh et al. 1998). 



5.2. Clusters of Galaxies 

The measurement of 6* g and r* g directly in the field is not possible given current data. However, the 
observed evolution of galaxy properties in clusters does give us a handle on the evolution of the distribution 
of galaxy formation. Here, we investigate cluster properties in the simulations. First, we describe the effect 
of the clustering of galaxy formation on cluster mass-to-light ratios, and how that subsequently affects 
estimates of ft based on mass-to-light ratios. Second, we discuss the Butcher-Oemler (1978, 1984) effect. 
We select the three objects in the simulation with the largest velocity dispersions, whose properties are 
listed in Tables |l| and || at redshifts z — 0.5 and z = 0. This number of clusters of such masses in the 



volume of the simulations is roughly consistent with observations (Bahcall 1998). Clusters #2 and #3 are 
about 5 h^ 1 Mpc away from each other at z = 0, and have a relative velocity of 1000 km s _1 almost directly 
towards each other; thus, they will merge in the next 5 Gyrs. 

To calculate mass-to-light ratios, we must calculate the luminosity associated with each galaxy particle. 
There are two steps: first, we assign a model for the star formation history to each particle; then, given 
that star formation history, we use spectral synthesis results provided by Bruzual & Chariot (1993, 1998) 
to follow the evolution of their rest-frame luminosities in B, V, and K. The simplest model for the star 
formation history of a single particle is that its entire mass forms stars immediately, the instantaneous 
burst (IB) approximation. However, it is unrealistic to expect that all of the star formation will occur 
immediately; at minimum, it will at least a dynamical time for gas to collapse from the size of a resolution 
element to the size of a galaxy. Thus, we also consider an alternative model based on the work of Eggen, 
Lynden-Bell, & Sandage (1962; hereafter ELS), which spreads the star formation out over the local 
dynamical time. We assign a star formation rate of 

dm* m n t 



-e 



dt ^dyn ^dyn 



"t/tdy 



(9) 



to each particle, where m p is the mass of the galaxy particle and idyn is the local dynamical time at the 
location of the particle at the time it was formed. The evolution of the mass-to-light ratios of clusters are 
qualitatively the same for both the IB and ELS models; the evolution of the B — V colors of cluster galaxies 
does differ somewhat between the models, as described below. 

We can calculate the mass-to-light ratios of these clusters in the B-band, Tg = (M / Lb) / (Mq / Lq,b) , 
and compare them to the critical mass-to-light ratio T b, c necessary to close the universe. We do the same 
in the V and if-bands, and show the results in Tables [l] and || at redshifts z = and z — 0.5, respectively. 
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An observer in the simulated universe would thus estimate Q cst (z = 0) = Tb/Tb jC ~ 0.4-0.5 based on 
cluster mass-to-light ratios, not far from the correct value for the simulation (Qq — 0.37). The galaxy bias, 
which causes fio.cst to underestimate f^o by a factor of about 1.6, is partially cancelled by the relative fading 
of the light of the older cluster galaxies with respect to the younger field galaxies. The bias effect (as we 
have shown above) increases in importance with redshift, while the differential fading of cluster galaxies 
decreases in importance. Thus, in the simulations, the cluster mass-to- light ratio estimates of f2 es t(z) 
decrease with redshift, even though O(z) itself increases. We find similar results for analyses performed in 
V and K, which we also list in Tables [l] and ||. For the ELS star formation model, the mass-to- light ratios 
are all higher, but the estimated values for f2 and their redshift dependence do not change substantially. 

Carlberg et al. (1996) studied 16 X-ray selected clusters between z ~ 0.2 and 0.5, in the K-corrected 
r-band. Their estimate of tt(z sa 0.3) = 0.29 ± 0.06 is not far from what observers in our simulation would 
conclude for that redshift, Sl(z ~ 0.3) « 0.38, using T^-band luminosities. However, the correct value in the 



simulations is fl(z = 0.3) = 0.56. Carlberg et al. (1996) found no statistically significant dependence of 



either T r or T,.. c on redshift, but their errors are large enough to be consistent with the level of evolution 
of these quantities seen in the simulations. We caution, therefore, that estimates of fio from cluster 
mass-to-light ratios can be affected by bias, by the differential fading between the cluster galaxies and the 
field galaxies, and by the evolution of both these quantities with redshift. 

As an aside, Tables [3] and || also show that Qq estimated from the baryonic mass in each cluster, 
QbM tot /Mb aTyons , is biased slightly high. This effect occurs because of the slight antibias (5 - 10%) of 
baryons with respect to mass at low redshift (due, presumably, to shocking and outflows in the gas which 
prevent baryons from flowing into the potential wells as efficiently as dark matter does). This result is 



consistent with previous work ( Evrard 199C ; Cen fe Ostriker 1993 ; White et al. 1993 ; Lubin et al. 1996 ) . 
Measuring Mbaryons is possible in the clusters because the gas is hot enough to be visible in X-rays; the main 
difficulty is in measuring f2b, which is generally based on high-redshift deuterium abundance measurements 
(e.g., Buries fe Tytlcr 1998j ) and Big Bang nucleosynthesis calculations (e.g., [Schramm fc Turner 1998 ). The 



fact that the baryon mass in clusters is closely related to the total mass makes this method of estimating 
fio less prone to the complications of hydrodynamics, galaxy formation, and stellar evolution than using 
mass-to-light ratios. 

From the results of Section ||, we expect that cluster galaxies at z = 0.5 should be younger and bluer 
than cluster galaxies at z = 0. Does this account for the Butcher-Oemler effect — the existence of a blue 
tail in the distribution of B — V colors of cluster galaxies at z ~ 0.5 (Butcher & Oemler 1978, 1984)? 
The simplest approach to this question is to ask whether the mass ratio of recently formed galaxies to all 
galaxies in the cluster changes from z = to z = 1. Figure [ll] shows the distribution of formation times 
of galaxy particles in each cluster at z = and at z = 0.5, assuming the IB model of star formation. For 
comparison, the thick histogram shows the star formation history of the whole box. At both z = and at 
z = 0.5, there has been no star formation in the previous 1 Gyr in any cluster. On the other hand, the more 
realistic ELS model allows star formation to persist somewhat longer. Figure [l2| illustrates this effect by 
showing the fraction of the stars in the clusters formed in the previous 0.5, 1. and 2 Gyrs at z = 0, z = 0.5, 
and z = 1, now using the ELS model. For comparison, we show the same curves for the whole box. Clearly, 
the fraction of young stars increases with redshift more quickly in the clusters than in the field. This is 
qualitatively in agreement with the observed Butcher-Oemler effect. 

We can look at the problem in a more observational context by calculating the B — V colors of cluster 
galaxies as a function of redshift. Following Butcher & Oemler (1978, 1984), we define the quantity /& (the 
"blue fraction" of galaxies) to be the fraction of mass in galaxy particles whose colors are more than 0.2 
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mag bluer than the peak in the B — V distribution. Observationally, /& ~ 0.05 at z = 0, but at redshift 
z = 0.5, f b ~ 0.1-0.25 (Butcher & Oemler 1978, 1984; |Oemlcr, Dressier, fc Butcher 1997|) . We consider both 



the IB and ELS approximations. Under the IB approximation, the colors redden as the galaxy populations 
evolve (about 0.1 mag between z = 0.5 and z = 0), but the shape of the B — V distribution changes little; 
for all the clusters at both z = and z = 0.5, fb ~ 0.05 for the IB approximation. However, using the 
ELS model improves the situation considerably; we find fb ~ 0.05 at z — and fb ~ 0.1-0.15 at z = 0.5, 
in good agreement with observations. Thus, these simulations may be revealing the mechanism behind the 
Butcher-Ocmlcr effect. 

It is worth noting that the astrophysics involved in the star formation history of cluster galaxies is 
clearly much more complicated than that which we model here. First, the star formation history inside 
each galaxy particle we create is surely more complicated than any the models discussed in the beginning 
of this section. Second, star-burst mechanisms such as ram-pressure induced star formation, which would 



affect galaxies falling into clusters ( Dressier fc Gunn 1990 ), or merger-induced star formation, which would 
affect galaxies in dense environments, could be important. These mechanisms are not modelled by these 



simulations. Third, as stressed by Kauffmann (1995), the population of clusters observed at z — 0.5 in an 
Q = 1 universe necessarily contains rarer peaks than the population of clusters observed at z = (Andreon 



fe Ettori 1999), and the differences in the formation histories of such different objects may account for the 
Butcher-Oemler effect. On the other hand, this effect should be unimportant for a low-density universe, in 
which structure forms relatively early. 



6. Conclusions 

We have examined the spatial distribution of galaxy formation as a function of redshift in a cosmological 
hydrodynamic simulation of a ACDM universe with Slo = 0.37. We quantified the evolution of the clustering 
of galaxy formation with the bias b*{z) and the correlation coefficient r»(z) of galaxy formation. In our 
simulations, r*(z) evolves considerably, whereas in the peaks-biasing formalism it remains unity at all 
redshifts. The time history of these quantities, combined with the gravitational debiasing, causes the 
evolution of b g (z) and r g (z) of all galaxies. As described in Section |^, this history of galaxy formation could 
be probed observationally by examining the evolution of &* g r* g , the linear regression of the star formation 
weighted galaxy density fields on the number-weighted galaxy density field. The mass-to-light ratios of 
clusters are also profoundly affected by this history, with the effect that while Vl(z) increases with redshift, 
the estimated 57 cs t(z) actually decreases. Finally, the evolution of the distribution of galaxy formation in 
the simulation bears a qualitative relation to the Butcher-Oemler (1978, 1984) effect, and under certain 
assumptions about the star formation history of each galaxy particle can even quantitatively account for it. 

Understanding these results completely will require future work. First, the low resolution of these 
hydrodynamic simulations (relative to state-of-the-art A-body calculations) calls into question the accuracy 
with which we can follow the process of galaxy formation, and makes us unable to identify individual 
galaxies in high-density regions. Thus we cannot follow adequately the effects of the merging and destruction 
of halos, especially in clusters. This fact was the basis of our decision to examine the stellar mass density 
field rather than the galaxy number density field. Observationally, it is possible to trace the stellar mass 
density crudely by weighting the galaxy density field by luminosity, especially at longer wavelengths. The 
evolution of bias in the luminosity-weighted galaxy density field is dominated by the three mechanisms 
described in Section ^. The evolution of bias in the nwm&er-weighted galaxy density field is additionally 
affected by the process of merging in dense regions. We do not have the resolution to treat this effect here; 
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Kravtsov & Klypin (1998) explore this using a high-resolution iV-body code. 



Equally important is to explore the dependence of the evolution of 6*(z) and r*(z) on cosmology. As we 
saw in Section their evolution was dominated by the evolution of the density-temperature relationship; 
as this relationship and its time-dependence varies from cosmology to cosmology, we expect 6*(z) and r*(z) 
to vary as well. Just as important is the uncertainty in the star formation model used, whose properties 
are evidently rather important in determining the clustering and correlation properties of the galaxies. 
Furthermore, understanding the statistical properties of galaxy clusters in the simulations will require more 
realizations than the single one we present. Finally, larger volumes than we simulate will be necessary to 
compare our results to the large redshift surveys, such as the SDSS, currently in progress. 

Since these simulations are expensive and time-consuming to run, we would like to find less expensive 
ways of increasing the volume of parameter space and real space which we investigate. It may be 
possible to capture the important aspects of the hydrodynamic code by using much less expensive iV-body 
simulations. One way to do so is to determine the relationship between galaxy density, mass density, 
and velocity dispersion at z = 0, and simply apply that relation to the output of an iV-body code at 



redshift zero, as suggested by Blanton et al. (1999). However, the final relationship between galaxies, 
mass, and temperature will depend both on the time evolution of the mass-temperature relationship, and 
the gravitational debiasing. For instance, if all galaxies form at high redshift, galaxies will have had time 
to debias completely. To take into account such effects explicitly, we want to apply a galaxy formation 
criterion at each time step, and follow the gravitational evolution of the resulting galaxies to find the final 
distribution at z — 0. If one could do so, one could run large volume dark matter simulations, in which 
one formed galaxies in much the same way as they form here. Motivated by Figure ^, we propose applying 
our measured (6*\6,T) on 1 ft, -1 Mpc scales in order to form galaxies in such an iV-body simulation, using 
the local dark-matter velocity dispersion or the local potential as a proxy for temperature. This approach 
assumes only that the dependence of galaxy formation on local density and local temperature is not a 
strong function of cosmology. This method would be similar in spirit to, but simpler than, the semi-analytic 
models which have been implemented in the past. 

The results of Section || suggest a simple way of exploring the effects of changing the galaxy formation 
criteria in such iV-body simulations, or for that matter in the hydro dynamical simulations themselves, 
as long as one is content with studying second moments and with ignoring the effects of the feedback of 
star formation. Given more output times than we have available from these simulations, one can simply 
determine the location of galaxy formation at each time step - after the fact - and use Equation (Q) to 
propagate the results to z = 0, or to whichever redshift one wants. Thus, one could easily examine the 
effects that varying the galaxy formation criteria would have. Compared to rerunning the full simulation, 
it is computationally cheap to determine the properties of (5* at each redshift (given the galaxy formation 
model) and then to reconstruct the evolution of b g (z) and r g (z). 

The continuity equation may also prove useful in comparing observations to models. For instance, 
large, deep, angular surveys in multiple bands, such as the SDSS, will permit the cross-correlation of 
different galaxy types as a function of redshift, using photometry to estimate both redshift and galaxy type 



( Connolly, Szalay, fc Brunner 1998 ). If one properly accounts for the evolution of stellar colors, one will be 
able to test the hypothesis that two galaxy types both obey the continuity equation. The sense and degree 
of any discrepancy will shed light on whether one type of galaxy merges to form another, or whether new 
galaxies of a specific type are forming. 



Estimates of the star formation rate as a function of redshift are now possible ( Madau et al. 199Q . 
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Pcrhaps the next interesting observational goal is to study the evolution of the clustering of star formation. 
The Butcher-Oemler effect tells us something about the evolution of star formation in dense regions, but it 
is necessary to understand the corresponding evolution in the field, as well. The correlation of star-forming 
galaxies at different redshifts with the full galaxy population may be a more sensitive discriminant between 
galaxy formation models than measures of the clustering of all galaxies. 
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A. Understanding Galaxy Formation in the Simulation 



The conditions required for galaxy formation in each cell of the simulation are: 



S > 5.5, 
m b > mj ee G" 3 / 2 Pb 1/2 C 3 

and 



S d + 1 Oi 
S b + 1 fib 



-3/2 



^cool ^ ^dyn 

V-v < 0. 



3tt 



32Gp' 



(Al) 



The subscript b refers to the baryonic matter; the subscript d refers to the collisionless dark matter. In 
words, these criteria say that the overdensity must be reasonably high, the mass in the cell must be greater 
than the Jeans mass, the gas must be cooling faster than the local dynamical time, and the flow must be 
converging. Assuming that the time scale for collapse is the dynamical time, we transfer mass from the gas 
to collisionless particles at the rate: 

Am* mb 
At ~ ~' 

where t* = t dyn if t dyn > 100 Myrs and t» = 100 Myrs if t dyn > 100 Myrs. 



(A2) 



To understand the effects of these criteria, consider Figure |T3|, which shows the distribution of 1 + 5 and 
T at the scale of a single cell at each redshift. We show the Jeans mass criterion as the diagonal grey line 
(points to the left are Jeans unstable), and the cooling criterion (for one-percent solar metallicity gas) as 
the curved grey lines (points in between the lines can cool efficiently). For these purposes we have assumed 
fib = S, 7 = 5/3, and /j, = tuh- Note that at high densities, satisfying the Jeans criterion automatically 
satisfies the cooling criterion. At high redshifts, there is a fair amount of gas which appears to have cooled 
into this regime, at 10 4 -10 5 K and at relatively high densities. As the density fields evolve, the densest 
regions become hotter, making it more difficult for the gas there to cool enough to allow galaxy formation. 
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Fig. 1. — Volume- weighted joint probability distribution of 1 + S, the mass density, and T, the gas 
temperature, for each output redshift. Contours are spaced logarithmically. Each field has been smoothed 
with a 1 h^ 1 Mpc top hat. Note that the range of possible temperatures at a fixed density increases with 
cosmic time. 
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Fig. 2. — The average rate of galaxy formation, 1 + 5*, as a function of 1 + 5 and T, for each output redshift. 
Each field has been smoothed with a 1 Mpc top hat. We measure the local rate of galaxy formation 
by examining the density distribution of "recently-formed" galaxy particles; that is, those formed in the 
previous 0.5 Gyrs. Contours are spaced logarithmically. Note that at any given density, the star formation 
declines with increasing temperature. As the densest regions become too hot to cool and collapse, they no 
longer become preferred sites for galaxy formation, which move to the relatively less dense regions of the 
universe. 
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Fig. 3. — Conditional probability P(l + + <5) at each of four redshifts. Greyscale stretch is logarithmic. 
Each field has been smoothed with a 1 hr 1 Mpc top hat. Solid line is the conditional mean (1 + <5*|1 + 6); 
dashed lines are la limits around that mean. Note that at later times there is increasing stochasticity and 



a turnover in (1 + 5«\l + S) at high density. While in Blanton et al. (1999) we showed where galaxies of 
different ages ended up at z = 0, here we show instead where galaxies form at each redshift. 
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Fig. 4. — The redshift evolution of the bias 6* (z) = cr* /<r (top) and the correlation coefficient r»(z) = 
{88*) /era* (bottom) "at birth," for top hat smoothing scales of 1 Mpc and 8 hr 1 Mpc. There is a strong 
trend with redshift. We also show, as the dotted lines, the peaks-bias prediction for 6*(z) and r*(z) on large 
scales, for M > 10 12 M Q . Note that r*(z) does not evolve at all in the peaks-biasing model. 
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Fig. 5. — Evolution of a burst of galaxies at z — 1, smoothed on 30 h^ 1 Mpc scales. The solid lines show 
the conditional mean galaxy density given the mass density, at z = 1, z = 0. 5 and z = 0. The dashed lines 
show the predictions at z = 0.5 and z = given the results at z = 1, according to the rule that A5 C = A<5. 




Fig. 6. — Same as Figure now smoothed on 1 ft, 1 Mpc scales and plotted on a logarithmic scale. 
Predictions are good in the regime: 3 < 5 < 100. 



-24- 




Fig. 7. — Evolution of the bias b c (R) and correlation coefficient r c (R) of a burst of galaxies at z = 1. Solid 
and dashed lines have the same significance as in Figure |^. In this case, predictions use Equation @ . 
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Fig. 8. — A comparison of the debiasing due to gravitational evolution and due to the shift in the location 
of galaxy formation. The solid line shows the evolution of galaxy formation, smoothed with an 8 /i" 1 Mpc 
radius top hat, from Figure |J. The dotted lines emanating from each output time indicate how the bias 
properties of galaxies formed at that time evolve according to the continuity equation. The top panel shows 
that gravitational debiasing acts more slowly than the shift of galaxy formation out of the high density 
regions. This is why old galaxies at redshift zero are still more highly biased than young galaxies, despite 
having had more time to debias. The bottom panel shows how the decline in r* tends to be washed out by 
the debiasing. 
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Fig. 9. — Evolution of the bias b g (z), correlation coefficient r g [z), and the combination b g (z)r g (z) of all 
galaxies on 8 h^ 1 Mpc tophat scales. Solid line is measured from the simulations. Dotted line is the 
reconstruction according to Equation (Q). The reconstruction suffers from the poor time sampling of our 
output. The dashed lines are the corresponding quantities for the recently- formed galaxies at each redshift, 
as in Figure ||. 
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Fig. 10. — Bias b„ g {z) = ct*/ct 9 , correlation coefficient = (S*S g ) /<j*<j g , and the linear regression 

b^ g r^g of the galaxy formation density field on the density field of all galaxies as a function of redshift. We 
show results for 1 h^ 1 Mpc and 8 h~ x Mpc radius tophat smoothing. The dotted lines are the peaks-bias 
predictions for the same quantities on large scales, comparing the bias of the mass in recently formed halos 
> 1O 12 M to that of all halos > 10 12 M Q . 
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Fig. 11. — Star- formation history in clusters at z = (thin solid line) and z = 0.5 (thin dotted line). 
Note that most star formation is complete in these clusters long before the universe is 7.8 Gyrs old, which 
corresponds to z = 0.5 (shown as the vertical dashed line). For this reason, there is no excess of blue galaxies 
at z = 0.5, although there may be at somewhat higher rcdshifts. For comparison, we show the star formation 
history of the full box as the thick solid line. This figure assume the instantaneous burst (IB) approximation. 
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Fig. 12. — Mass fraction of recently formed stars in the clusters (solid lines) and in the whole box (dashed 
lines), using the ELS star formation model described in the text. We show results for stars formed in the 
previous 0.5, 1, and 2 Gyrs. In each case, the rapid decline in the clusters compared to the held is evident. 
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Fig. 13. — Volume-weighted joint distribution of mass density and temperature at 0.2 hr 1 Mpc scales (the 
grid cell size of the simulation), at each output redshift. Greyscale and contours are as in Figure [I]. Note that 
at these scales, the peak in the cooling curve at 10 4 to 10 5 K is noticeable as a ridge at around 1 + 5 ~ 100. 
We overlay the Jeans criterion for galaxy formation as the diagonal line and the cooling criterion (for one- 
percent solar metallicity) as the two vertical curves. The simulation produces galaxies above the Jeans line 
and in between the two cooling lines. At all redshifts, the densest regions automatically satisfy the cooling 
criterion as long as they satisfy the Jeans criterion. 
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Table 1. Clusters at z = 0, at which time fl = 0.37. 





M tot (/i" 1 M ) 


(v 2 ) 1/2 (km/s) 


T B /T B , C 


Tv/T v ,c 


T K /T K , C 


flboMtot 1 Afbaryons 


#1 


9.847 x 10 14 


1076.505 


0.414 


0.366 


0.316 


0.390 


#2 


7.919 x 10 14 


951.378 


0.442 


0.401 


0.359 


0.401 


#3 


5.448 x 10 14 


832.657 


0.535 


0.498 


0.481 


0.409 



Table 2. Clusters at z = 0.5, at which time fl = 0.66. 





M tot (h- 1 M ) 


(v 2 ) 1/2 (km/s) 


T B /T B , C 


r v /T v , c 




^boMot/Mbaryons 


#1 


2.698 x 10 14 


751.243 


0.321 


0.275 


0.231 


0.716 


#2 


6.468 x 10 14 


1135.564 


0.451 


0.390 


0.321 


0.719 


#3 


3.839 x 10 14 


815.936 


0.397 


0.349 


0.304 


0.663 



